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LMIs decomposition approach. Next, the observer’s gain matrices are determined 
Lyapunov theory using the Lyapunov theory and the constraints are defined as linear matrix 
Takagi-Sugeno singular model inequalities (LMIs). Finally, a numerical simulation on an application 


example is given to demonstrate the usefulness and the good performance of 
the proposed dynamic system. 
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1. INTRODUCTION 

Over the last decades, the increase in the performance of equipment in terms of production quality 
and gain in productivity was accompanied by the complexity of the equipment. However, the presence of 
abnormal changes due to actuator or system or sensor faults can degrade system performances, hence the 
need to integrate fault detection and diagnosis (FDD) tools [1], [2] to maintain, for a long time, the desired 
performance of the whole system in various sectors. In particular, FDD has a very important role in 
monitoring the behavior of system variables and revealing faults, and it is performed based on the relative 
information to the system and its equipment. This information can be obtained by adding sensors to acquire 
measured states or observers to estimate non-measured states requiring expensive or difficult sensors to 
maintain. The state modeling of process dynamics is often obtained based on its state variables linked 
together by mathematical equations. If these processes have constraints, then it is necessary to use static 
equations to sufficiently characterize the studied process. Such systems composed of static and dynamic 
equations are called singular, or descriptor or implicit systems [3]. Recently, the FDD problem for singular 


Journal homepage: http://ijeecs.iaescore.com 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 Oo 173 


systems has attracted much attention in various fields such as mechanical engineering, computer science, 
civil engineering, electrical engineering and automation. 

Various techniques for detecting and estimating faults have been proposed for the class of linear 
systems [4], [5], and for the class of nonlinear systems [6]-[8] allowing to provide a closer representation to 
the real system, but which are difficult to exploit. Due to this complexity, it has become essential to work 
with a precise class of nonlinear systems such as Lipschitz systems, uncertain systems, bilinear systems or 
others. From these classes of nonlinear systems, we find the class of Takagi-Sugeno (T-S) [9] nonlinear 
systems, in ordinary or singular form. It has been introduced to compromise between the good precision of 
the nonlinear behavior of the studied system, and the use of techniques adapted to linear systems due to the 
convex sum property of its activation functions [3], [10]. 

There have been many methods of FDD [2], [11], [12] which can be classified into signal-based 
approaches [13], knowledge-based approaches [12], [14], and process model-based approaches [12], [13] 
which contain state observer-based method representing an analytical method having achieved several results 
in this field, and which depends on the mathematical model of the studied system without needing other 
components. Many publications have been interested in the design of observers for FDD [15]-[24] and have 
presented fruitful results. A residual generator for detecting and isolating actuator faults for a class of T-S 
fuzzy bilinear system is developed in [18]. Developing a novel fuzzy FD observer for FD of sensors faults of 
T-S fuzzy systems is the aim of the work presented in [19]. In [20], depicted a T-S unknown input observer 
to simultaneously estimate the interval of states and actuator faults for a class of T-S explicit systems. 
Another technique based on a robust fault estimation observer has been introduced for estimating actuator 
faults for a class of discrete-time singular systems [21]. In [22], a design of an adaptive observer is proposed 
for detecting sensor faults of an industrial servo system. For the fault diagnosis and reconstruction of the 
faults affecting the states of the system, in [23] suggested a new augmented linear parameter-varying (LPV) 
observer for a class of LPV models. The design of a combination of reduced-order LPV and full-order LPV 
unknown input observers, respectively, for FDD of actuator and sensor faults of industrial processes is 
presented in [24]. 

Most of these observers are synthesized to estimate only actuator or sensor faults while guaranteeing 
asymptotic convergence for various class of nonlinear system in continuous or discrete-time. The goal of our 
work is not to compare our approach with those already carried out, but rather to extend our results from the 
case of singular linear models [25] and T-S singular models with measurable premise variables [26] to the 
case of T-S singular models with unmeasurable premise variables while ensuring an exponential 
convergence, and simultaneously estimating the unmeasurable states and the faults at the level of actuators 
and sensors. In this work, for simultaneous estimation of states and faults, the novel suggested technique 
consists to associate for each local model a local observer. Then, the proposed fuzzy observer is obtained by 
an aggregation of the local observers. Our contribution is based on the separation of the dynamic equations 
from the static equations which makes it possible to facilitate and minimize the computation by obtaining the 
static states just from the dynamic states already found. The design conditions are expressed in terms of 
LMIs. This observer is applied for both actuators and sensors faults for a class of T-S singular model in the 
case of unmeasurable premise variables. The paper is composed of five parts that are presented as follows: 
Section 2 exposes the class of the studied system. Section 3 provides the synthesis of the proposed observer 
and the stability conditions. The numerical results of the application example are given in section 4. Section 5 
is devoted to a brief conclusion. 


2. MATHEMATICAL FORMULATION OF THE CONSIDERED MODEL 
In this paper, the following class of continuous-time Takagi-Sugeno singular model (CTSSM) with 
unmeasurable premise variables in presence of actuator and sensor fault is considered (1), 
f Mż = Xi, pi(B) (Aiz + Bit + Maiva) (1) 

eS D pil(B) (Ciz + Dit T DaiYa + MsiVs) 


where z = [Z] Z7]! e€ R” is the state vector with Z4 € R” is the vector of dynamic variables, Z, € R"" is 
the vector of static variables, t E R? is the control input, y € R” is the measured output vector, vg € R° 
and v, € R? are the actuator and sensor fault vectors, respectively. The matrices A;, Bi, Mais Ci Di Dai and 
M,; are real known constant matrices with adequate dimensions related with the it” local model with, 


_ (Ani Arni ; — (Bi). — (Mai). — 
a= (1E 22) B= (G) Ma = (KE): C= Cu Co Q) 


2 
Mai 


State and fault estimation based on fuzzy observer for a class of Takagi-Sugeno singular... (Kaoutar Ouarid) 


174 m) ISSN: 2502-4752 


the rank of the matrices Az ; are equal to n — r and it is supposed to be invertible. q represents the number of 
sub-models, and the premise variable f is supposed to be real-time accessible. 


{ Mz = A;z + Bit + Maida 3 
y= Ciz + Dit + DaiVa + M;iVs ( ) 


The transition between the contributions of each sub model (3) is ensured by the terms p;(8) which 
represent the weighting functions, depending on the states of the system and verifying the property of the 
convex sum, 


ELAB) =1;0<p:()<1;i=1,..,q (4) 


the matrix M whose rank (M) = r < n is assumed to have the following form, 


"-( 9) o 


Assumption 1: Assume that [3]: 
— (M, A;) are regular, i.e. det(sM — A;) # 0 VseC 
— The sub-models (3) are impulse observable and detectable 

The separation of the dynamic equations from the static equation in each sub-model (3) is the aim of 
our approach, and then the aggregation of the resulting sub-models allows obtaining the global fuzzy model. 
So, using the expression of the matrices (2) and (5), the sub-model (3) can be written in the following second 
equivalent form [3], 


Żı = AyyiZ, + Ay2iZ2 + Byit + Maia 
0 = A2121 + A2222 + Bit + Miia (6) 
y= CiiZa + CziZ2 + Dit + DaiYVa + M;iVs 


by finding the expression of the static variable Zz, and replacing it in (6), we obtain, 
Ži =JiZ, + Lit + Naiva 


Zz = SiZ, + Vit + Raia (7) 
y= PZ, + T;T + KgiVa + MiV; 


where, 
Ji = Ari + A21 
Li = Bii + A12:;V; 
Nai = Mai + A12iRai 
Si = Azz Aati 
V; = —Az7;Boi (8) 
Rai = —Az3iMii 
Pi = Cii + CaS; 
T; = Di + CaiV; 
Kai = Dai + CziRai 
let define, 


v=(%) ©) 
which is equivalent to the following state representation, 
Z, =JiZ, + Lit + Nv 


Zo = SZ, +V;t + Riv (10) 
y = PiZi +Ti;t + Kiv 
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N; = (Nai 0) 


R; = (Rai 0) (11) 
Ki = (Kai Msi) 


then, from (10) p;(B) can be rewritten as, 
Pi(B) = pi(Z1, Z2 = S:Z1 + Vit + Riv) = pi (Z1 T, v) = pi (w) (12) 
with w? = [Z] t" v]. So, the system (1) can be rewritten under the following equivalent form, 


Z1 = DL, pi(w) Q:Z1 + Lit + Niv) 
Z = DE, pi(w) (SZ, + Vit + Riv) (13) 
yY = LL, pi() (PZ, + T;t + Kv) 


Assumption 2: Assume that v is considered in the following form, 
V = Co + Cat + Cyt? +o + Cm tS (14) 


where cp ; k = 0,1, vey Mp are real unknown constant parameters and the (my + 1) time derivative of the 
fault is null. 


Let, 
Pk =v with k = 1, ... mp +1 (15) 
then, 
o ey. with k = 1,...,m, (16) 
mf 


thus, we rewrite the system (13) under the equivalent augmented state form as follows, 
X= DE, pi(9) QX + Lit) 


X = LE, pi(8) (iX + Viz) (17) 
y = DE, pi(O) (PX, + 7:7) 


where, 
XT =(ZI Pi o Pmp) 
Xz = Z 
~ (% 
a) 
Ji Ni 0 > 0 
0 0 FJ + 0 
Jes p a o n (18) 
0 0 0 > I 
0 0 - 0 0 
brEt 0 0 oy 
Š =(S; Ri 0 0) 
P =(P; K; 0 0) 


3. RESEARCH METHOD 


The following section shows the design of new structure of fuzzy observer allowing the 
simultaneous estimation of the unmeasurable states and unknown faults of the equivalent structure (17) of the 
CTSSM (1), 
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X, = Di plô) GX, + Lit- Gi — y)) 
Xe = Lies Pi) (SiX, + Vit) (19) 
f = Xi pi) (PX, + Tt) 
such that the estimated vectors of (X4, X2) and y are denoted by (X,, Xz) and Ẹ, respectively. 
For i = 1,---,q the term G; expresses the observer gain for the it” submodel such as the estimated 


of the augmented vector of the states and faults tends asymptotically towards the real vector. 
Defining, 


eo) k) any 


substituting (17) and (19) into (20) gives the following static and dynamic equations of the state estimation 
error, 


[a = D i) Ji + Lit — GO — y)) — Ei p:(0) UX, + Lt) (21) 
ez = YE pi) SX, + Vit) — DL, Pi) (iX + V7) 
equivalent to, 
k = YE, Pi) Jier — Gi — y)) — LL (0:9) — p: (®)) JX + Lit) (22) 
ez = Dk, pi(9) Siei — DIL, (i (0) — 9:(8)) (SX + Vit) 
let consider F; = J;, Li, Si, V; and, 
EL (i (0) — p:(6)) Fi = Dij pi(O)p;(8) A Fij (23) 
Then, by using the expression (23) the system (22) becomes, 
6, = YL, pi(8) Gie - Gi -— y)) — Èi j=1 Pi(9)p;(6) (Afi, X + AL;;T) (24) 
ez = YL, pi (Ô) Sie, — Èi j=1Pi (0) p;(8) (ASX, + AV;;T) 
as ys pi(0) = 1, we obtain, 
ė& = Xia p:(0)p;(®) J;e = GV -y))-— Dai pi(0)p;(6) (Af F AL;;T) (25) 
en = Èi j- Pi) 0; 6) (Še, = AS, X — AV;;T) 
in the same way, we can get y and ¥ as follows, 
P =P piCO)Pn (Ê) (CP, + APin)X1 + (Th + ATin)T) (26) 
P = Viner Pi(9) Pn (8) (PrX + ThT) 
with AP,, = P, — P, and AT;, = T; — Ty. By the substitution of (26) in (25), we get, 
é = Dipha pi (9) p;(8)pn(8) (ines + Fijn Xz + Vijnt) (27) 
e2 = ae pi()p; (8) (Sje, — AS;;X, — AV;;7) 
with, 
Zijn = GAPin — Ajij (28) 


Vijn = GATin = AL; 
ij,he (1,-,q) 
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therefore, to demonstrate the convergence of e towards zero, it suffices to demonstrate that e, converges to 
zero. Considering & = (e7 XT)", we get, 


{* Dae p:(8)p;(8)on(9) (Zijn 6, + Vint) (29) 
e1 = Që, 
with, 
mC 
rn (1) o 
Q=( 0) 


guaranteeing the stability of (29) while attenuationg the effect of Tt on e; is linked to the determination of the 
observer gains G; fori = 1,---,q. 


Theorem: Under assumptions 1 and 2, if for the CTSSM (1) there are matrices P,,P,,W; for 
i = 1,---,q, and a positive scalar € for a given ð > 0 which satisfy the LMIs (31), then it will be possible to 
determine the observer gains, that ensure the exponential convergence to zero of the estimation error, 


6; Yijn Pijn 


Yin 4i Poli | <O0V(ij,h) € (1,0)? GD 
ijh [i P, —é 


with, 


Sin = JPP, + Pij; — PR Wj? — W;P, + 20P, +1 
Yijn = WP; — Pa) — P Qi- J) 


ahs Se (32) 
Pijn = WT; — Th) — Pi (Li — Lj) 
the gains of the observer G;,i = 1,---,q are obtained by, 
Gi = PW, (33) 


the attenuation level is, 


Proof of Theorem: Let us consider the following quadratic Lyapunov function as follows, 
V(é,) = êP é,,P =P’ >0 (35) 


with, 
_ (Py 0 
P=(0 a) em) 
the derivative of V with respect to time is, 
V = Yo iner1PiC);()pn(8) (Ef (ZfinP + Pin) č + EL PVijnt + 07 Yi, P ê) (37) 


to guarantee the stability of (29) and the boundedness of the transfer from the input T to e,, 


led < @ ijti] #0 (38) 


IItll2 


we consider, 
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V+ele,—a’t™t<0 (39) 
so, the exponential convergence of the estimation error is guaranteed if, 
V +e? e — æ@?t" t < —20V ô > 0 (40) 


inserting (29) and (37) into (40) leads to the following inequality, 


V= Ei nas Pi (0)p;(@)pn(6) (Ef t) ijn (Es. T)" < 0 (41) 

with, 
ZinP + PZijn + Q7Q+20P PYijn 

lin = ( WIP e Ha 
then, the inequality (41) is contented if, 

Tijn < 0 Yi j,h € (1,+,q) (43) 
taking account (28), (30), (36), and the following change of variables, 

o “ 


we can deduce the LMIs (31) presented in the Theorem that complete the proof. 


4. RESULTS AND DISCUSSION 
To display the benefits of the suggested observer, we consider the following CTSSM which is 
affected by faults, at the level of actuator and sensor, and subjected to unmeasurable premise variable, 


= Di-1 pi(B) (Aiz + Bt + MaYa) (45) 


y = Cz + M,v, 


where z E€ R*,t E€ R,y E RŽ, va € Rand v, E€ R are the vectors of states, input, output, actuator fault 
and sensor fault, respectively. 


0 1 0 0 0 1 0 0 
_{-2.5 -0.75 0 0.025) , _[-—2.696 -0.75 0 0.025 
ala 0 1 —0.4 0 E 0 1 —0.4 0 (16) 
-2.5 -0.75 0 0.075 —2.696 -0.75 0 0.075 
0 
K 0 i ara Aa Ot Sk OW ana 
iz 0 ; Mg = B; C= (5 1 0 1) Ms = (6) a) 
—0.125 


pı (B) and p3 (B) represent the weighting functions, 


p(B) = n 
max min (48) 


p2(B) = ~fmax-F 


Bmax-Bmin 


where the expression of the premise variable £ is, 


-5-5z2 
p = = E [Bmin Bmax] (49) 
in order to apply the suggested fuzzy observer (19) on our application example (45), it suffices to represent it 


in its equivalent form (17). Thus, by using the Theorem with ð = 0.1, we obtain the following observer gains 
G and G3, 


Indonesian J Elec Eng & Comp Sci, Vol. 25, No. 1, January 2022: 172-182 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 Oo 179 


—16.2542 0.4992 —16.3254 0.5089 
60.4899 19.3705 59.6945 28.7403 
G = 443.5258 136.5402 |. G= 437.8298 202.7100 (50) 
= 77.1069 23.0732 |’ ? 76.1530 34.2843 
35.2472 —5.6930 35.5444 —7.7807 
12.2730  —3.6819 12.4475  -—5.1612 
the simulation results are given in Figures | to 4 where the input signal is given by, 
_ (t-—2whent <2 
= { 0 else On 


Under Assumption 2, the trajectories of actuator and sensor fault signals, which are applied respectively 
during the intervals [40, 160s] and [200, 320s], their first order derivatives, and their estimates are shown in 
Figures 3 and 4. 

These results demonstrate that the suggested fuzzy observer gives good performances in estimating 
unmeasurable states and unknown faults while catching up with unwanted variations. This approach has the 
benefit of being applied at the level of a large class of nonlinear systems. This is due to the fact that it is not 
required to know the value of the Lipschitz constant that can influence the resolution of LMIs [27], as well as 
without being limited by the condition of the rank between the matrices such as in [28]. 


0.25 T T T T T T 
— z: Real 


— 3z: Real 


xi Estimated ---2Z,: Estimated 


— 2, : Real 


==: Estimated 


Figure 2. Z3 and z4 with their estimates 
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— ~y: Real 


--—- V, : Estimated —--%, : Estimated 
a 


— V,: Real 
---V,: Estimated 


0 50 100 150 200 250 300 350 400 à 0 50 100 150 200 250 300 350 400 


Figure 4. v, and ú, with their estimates 


5. CONCLUSION 

This work is addressed to the design of fuzzy observer for simultaneous estimation of unmeasurable 
states and unknown faults, for Takagi-Sugeno singular models in continuous time. The main idea of this 
paper is to extend the results developed in the case of measurable premise variables. The diagnostic 
procedure is based on the separation of static equations from dynamic ones. Using this, the determination of 
the static variables will be deduced from the computation of the dynamic variables. At last, an example of 
application is presented in order to highlight and confirm the effectiveness of the proposed approach in the 
estimation of the states, and the faults of actuators and sensors. 
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